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§ 1. Introduction. 

The problem of the eqiiilibriiim of an elastic solid under given applied forces is one 
of great difficulty and one which has attracted the attention of most of the great 
applied Mathematicians since the time of Euler. Unlike the kindred problems of 
hydrodynamics and electrostatics, it seems to be a branch of mathematical physics in 
which knowledge comes by the patient accumulation of special solutions rather than 
by the establishment of great general propositions. Nevertheless, the many and 
varied applications of this subject to practical affairs make it very desirable that these 
special solutions should be investigated, not only because of their intrinsic importance 
but also for the light which they often throw on the general problem. One of the 
most powerful methods of the mathematical physicist in the face of recalcitrant 
differential equations is to simplify his problem by reducing it to two dimensions. 
This simplification can only imperfectly be reproduced in the Nature of our three- 
dimensional world, but, in default of more general methods, it provides an invaluable 
weapon. 

It was shown by Airy^ that in the two-dimensional case the stresses may be 
derived by partial differentiations from a single stress function, and it was shown 
latert that, in the absence of body forces, this stress function satisfies the linear 
partial differential equation of the fourth order V*x = 0, where V^ = V^ . V^, and V^ 
is the two-dimensional Laplacian d^/dx^ + d^fdy^. 

It might have been expected that these results would have opened the way for 
a theory of two-dimensional elasticity of the same generality as the two-dimensional 
potential theory. This has not, however, been the case. This is due in part to the. 
greater analytical difficulties which attend the discussion of the two-dimensional 

* ' Brit. Assoc. Rep.,* 1862, p. 82. 

t W. J. Ibbetson, * Proc. Lond. Math. Soc.,' vol. xvii., 1886, p. 296. For a history of this part of the 
subject see Love's * Elasticity,' 2nd edition, p. 17. 
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solutions of V*x = ^^ ^^^ compared with V\ = 0. The analogues of many of the 
important properties of the simpler equation have yet to be discovered if they exist at 
all. Some progress has been made, and in this connection we may mention the work 
of J. H. Michell'^' who estabhshed a general theory of inversion which, with some 
important differences, follows the potential theory fairly closely. 

No doubt the analytical difficulties have been the chief obstacle to progress, but 
perhaps the theory has not in recent years received attention which it would have 
received but for a certain physical difficulty. A truly two-dimensional elastic system 
is not so easy of realisation as might seem to be the case at first sight. If the stresses 
are everywhere parallel to the ayy plane and independent of z there will in general be 
a varying displacement parallel to z. If the displacements are everywhere parallel to 
the xy plane and independent of z this can only be secured by the application of 

a stress zz which varies from point to point and is perpendicular to the xy plane. This 
difficulty was in a large measure removed by a theorem established by Filon, which 
has been called the theorem of generalised plane stress, t It states that if the average 

value of the stress zz be taken throughout the thickness of a plate parallel to the xy 
plane, then the ordinary two-dimensional theory will give accurately the average stresses 
through the thickness of the plate if the elastic constants of the material are modified. 
If X, JUL denote the true elastic constants, X must be replaced by X^ = 2X/i/(X + 2/x) while 
JUL remains the same as before. This theorem attains an even greater importance when 
considered in the light of Miohkll's theorem,! that if a plate bounded by any 
number of bounding curves is in equilibrium under forces in its plane applied over the 
boundaries, then, provided the forces applied over each boundary taken separately 
are in equilibrium, the stresses are everywhere independent of the elastic constants. 

The hypothesis that the average value of zz vanishes throughout the plate, 
while certainly not accurately true in the majority of cases, will probably give 
a very close approximation in the case of a thiri plate where parallel faces are 
unstressed. 

In the light of this generalisation it is of considerable importance that the two- 
dimensional problem should be worked out more thoroughly. The two-dimensional 
solutions of \^\ = have been investigated in several systems of curvilinear co- 
ordinates. Owing to the special importance of the problem of the rectangular beam 
the solutions in Cartesian co-ordinates have naturally received a considerable amount 
of attention. Michell gave the general form of the stress-function in polar co- 
ordinates, thus opening the way for the solution of the problem of a plate bounded by 

* **The Inversion of Plane Stress," ' Proc. Lond. Math. Soc.,' 1901, vol. xxxiv., p. 134. Many of the 
results of the present paper can be obtained by an application of Michell's methods, but it has proved 
more convenient to proceed on different lines. 

t 'Roy. Soc. Phil Trans.,' A, 1903, vol. 201, pp. 63-155. 

I * Proc. Lond. Math. Soc./ vol. xxi., 1900, p. 100. 
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two concentric circles, or an infinite plate containing a circular hole under any- 
given tractions applied over its boundaries. In his lectures at University 
College, London, in 1912, Prof. Filon gave the complete solution of this problem 
determining the stresses and displacements when the stresses on the boundaries 
are expanded in Fourier series, and I am not aware that this solution has ever 
been published. An outline of the solution in elliptic co-ordinates is given in Love's 
' Elasticity/^ 

In this paper the complete solution is given for bipolar co-ordinates, for which the 
co-ordinate curves are co-axial circles. This solution enables us to treat the problems 
of an infinite plate containing two circular holes, a semi-infinite plate bounded by a 
straight edge and containing one circular hole, and a circular disc with an eccentric 
circular hole. 

In the second Section the equations are expressed in bipolar co-ordinates and 
formulae are established for the displacements in terms of the stress-function. 

In the third Section the stress-function is obtained in a convenient form and the 
terms giving rise to many valued displacements are separated out. 

The fourth Section is devoted to the determination of the coefficients in the stress- 
function when the tractions over the boundaries are given in Fourier series, and to an 
examination of the convergence of the resulting series. From the results established 
in this section it appears that the solution is complete, for the stress-function can 
always be uniquely determined when the tractions are given, provided that the 
applied forces taken as a whole are in equilibrium. 

The remaining sections are occupied with the examination of some of the simpler 
applications of the theory. Section 5 gives the solution for a circular disc with an 
eccentric hole (or a cylinder with eccentric bore) when the two boundaries are under 
different hydrostatic pressures. It is foimd that the solution of this problem can be 
expressed in finite terms. An important particular case of this problem is discussed 
in Section 6, namely, a semi-infinite plate with a straight, unstressed boundary and a 
circular hole under a uniform normal pressure. This will give the stresses near 
a rivet hole while the hot plastic rivet is being forced home under pressure. 
This solution is interesting from another point of view, for if the ratio of the 
radius of the • hole to its distance from the edge is suitably adjusted, the point of 
greatest tension will be on the straight edge while the point of greatest stress 
difference is on the circular boundary. It thus suggests a crucial test for the 
rival theories of rupture, — the greatest tension theory and the greatest stress- 
difference theory. 

Section 7 deals with a semi-infinite plate with an unstressed circular hole 
under tension parallel to its straight edge. The solutions are in the form of infinite 
series, but the more important aspects of the problem are illustrated by numerical 
tables. 

* 2nd edition, p. 259, 
2 P 2 



268 



mi, G. B. JEFFERY ON PLANE STRESS AlSfft 



§2. The Co-ordinates. 
Let us take curvilinear co-ordinates defined by the conjugate functions 

a •\-i8 = loQ' rW L .........(]) 

where x, y are Cartesian co-ordinates and a is a positive real length. Solving for 

x^ y^ we have 

a sin /3 a sinh a 



X ~ 



cosh a— cos B 



y 



cosh a— cos /5 



»« • • « « • A <w / 



Elements of ai^c measured along the normals to the curves a, ^ - 

respectively c!a//i, (!^/A, where 



constant are 




a 



a 



\da, 



from which we have 



h = (cosh a— cos ^)/a, . 



(3) 



The general scheme of co-ordinates is shown 
in fig. 1. If Oi, O2 are the points 0, —a and 
0, a respectively and P any point in the plane, 
and if the radii from Oj, 0^ to P are of lengths 
Ti, 7*2 and are inclined at angles 6^, $2 to the 
axis of X, then a = log r^/r^ and /3 = di—O^. 
The curves a = constant are a set of co-axial 
circles having Oi, Og for limiting points. "The 
circles corresponding to positive values of a lie 
above the a^-axis and those corresponding to negative values below, while the i:c-axis 
itself, which is the common radical axis, is given by ct = 0. The curves /3 = constant 
are circles, or rather arcs of circles passing through Oi, O2 and cutting the first set of 
circles orthogonally. On the right-hand side of the t/-axis /3 is positive and on the 
left-hand side negative, while on the |/-axis 8—0, except on the segment OjOjj, 
where /3 = ±7r. At infinity a = 0, ^ = 0, and at Oj, O2 we have a = — 00 and + 00 
respectively. 

We have thus a set of co-ordinates adapted for the consideration of two-dimensional 
problems in which the region considered is— 

( 1 ) A finite region bounded internally by a circle and externally by a larger and 

non-concentric circle. 

(2) A semi- infinite region bounded externally by a straight line and containing a 

circular hole. 

(3) An infinite region containing two circular holes of any radii and centre 

distance. 
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If the displacements in the directions normal to the curves a and ^ constant are 
u, V respectively, the strains are given bj^*^ 



^aa — f^ 



d-u 



V 



ah 



da 3/3 



^^^/3 



h 



dv 



oh 
a a 



a 



a 



ca 



(hv) + 777; (h/ii), 



^^ 



and the corresponding components of stress by 



m 

a/3 



'-= \{eaa + e^p)+2jLie 



act? 



X(6„«4-6^^)+2/x%, 



• t 



M^a^• 



j 



• • • • » » » X'x. J 



These stresses may be derived from a stress-function, so that in rectangular 
co-ordinates 






3= 



X 



dy 



2' 



ax ay 



yy 



9^x 



Transforming these equations to curvilinear co-ordinates we obtain 



aa = 



A^U^^ 







h 



dhd 



X 



"\ 



3a 3a 



/3/3 = /^ 



a/3 



3 



9 6, be 



a 



// 



-^2 






(hx) 



hikk, y 



3/3 3/3 



+ h 



dVi 



X 



3a 8,8 '^ 3a 3/3" 



• • • » t t » 



-/ 



(5) 



We will usually find it convenient to deal with hx instead of x itself, and in our 
particular co-ordinates these equations become 



aaa = 



r ^2 

< (cosh a— cos /3) 7^7^^ 



sinh 



a 

Oa 



sin ^ — + cosh a [ (Ax)- 

a 



^ 



3^ 



a/3/3 = -^ (cosh a— cos /3) 7^— ; — sinh 



V. 



^ 2 



a 



o 
da 






sm/3™--h cosh^j>(/ix)3 ^ . . (6) 



aa^ = — (cosh a— cos ^) 



da 8/3 



J 



We may note that 



a (act— /3^) = (cosh a— cos ^) 



^2 



32 



.3/8^ aa 



~2 +1 (%)j ..... (7) 



-)f 



Love, * Theory of Elasticity/ p. 54. 
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so that if hx and its second differential coefficients are finite at infinity (a = 0, /3 = O) 

we have there aaa = a/3/3 = /?x and aj8 = 0. 

In the absence of body forces tlie stress-finiction satisfies V\ = 0. In curviHnear 
co-ordinates we have 



and, taking hx as the dependent variable, we have in onr co-ordinates 

aV\ = \ (cosh a— cos /3) (tt^ + ;^) -2 sinh a ;^ -2 sin B :^+ cosh a+ cos ^K^^x)• 

I \da^ 0/37 Ba d/3 j v /v/ 

Eepeating the operator, a little reduction leads ns to the following transformation 
for V^x = • 

/ rl^ rl^ r)^ r)^ '^^ \ 

Thus by considering hx instead of x we have a Unear equation with constant 
coefficients. 

Before proceeding to the discussion of its solutions, we must investigate the method 
of determining the displacements corresponding to a given stress- function, in order 
that we niay ascertain whether and under what conditions these are single-valued. 
This is particularly necessary in our case, as one of the co-ordinates, ^, is itself many- 
valued. 

Adding and subtracting the first two equations (4), and leaving the third as it 
stands, substituting for the stresses in terms of the stress- function, and for the strains 
in terms of the displacements, we obtain the following three equations : — 



|.|,= | + V»}-|{*'|+2M.} = 0. (10) 

From the last two of these it appears that we may define a new function P such that 



ap 



^ =^^|~ + 2ya^V, ......... (12) 

VOL dp 

r~r ^^^^^ /?' 7^ -p ZifXflU^ . , , , . , , , , liO/ 

dp ca 

V^P = 0, 
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and we have still to satisfy (9). Substituting for it, v in terms of P we have 



h 



2 O ( 

cm 



1 aP\ . ,,2 8/1 3P\ X + 2fi y2 



:h'?:l3j'^^'''d^\K'da 



\ + /UL 



X'^ 



which may be re-arranged thus 



calx + M Sot A^o/3j 3/3 1 X-fM 3/3 A'^ 3a 



It follows that a function Q exists such that 



h 



h 



aa 



a 



X + /X 3/3 '■ 






d 



5P >H:_2^^ yj2 "x 

9/S X + /X ra 



0. 



(14) 



(15) 



Eliminating P by differentiating with regard to /3 and a respectively, and adding, 
we have 



cr {hQ) ,, _d%__ __ X4-2M II 3_ U2 8x\ 1 _i / p axV 



8ct 3/3 



oa c 



/3 ' 2(X + m) l^rJaV" a J hd^\ ' 3/3/ 



which becomes in our co-ordinates 



d' (kQ) 






aa3/9 2(X + m)\ 3«' 



3/3^ 



• f> 



. . (16) 



There is, however, a further condition to be satisfied by Q corresponding to the 
condition VT = 0. Differentiating (14) and (l5) with regard to a, ^ respectively, 
and subtracting we have 



a ' 



oQ\ 






8 /j,^ cQ '^ 



X 4- 2/x d 



'^ 



d^\ dpi X + M 13«\ 3,8/ 3,8 






or in our co-ordinates 



3^ 



3^ 



o« op 



2(X + 2M)3-'(/iQ) 

X-f/x 



8ot 3^ 



. . , . (17) 



These two equations connecting Q and x ^re consistent, for, if we eliminate Q by 
appropriate differential operators, we have 



c 



3^ 



da' d/3 






which is readily seen to be identical with the condition V\ = ^' ^^ given in (8). It 
is obvious that hQ, satisfies the same differential equation, and hence it also is a 
solution of V^Q = 0. 
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We haye therefore from (16) 



Z(\-^a)jJ{ COL vfT J 



2 (XH-/i) J J t oa^ 0/3 

and from (12), (13) and (14), (15) 



^ 



3Q 



2^«==^/,,£X„/,£^^ . . (19) 



OY . 1 gQ 



2Mt' = --^- // ^ + A P (20) 

X-{-fx op aa 

It is readily seen that these equations determine u and v apart possibly from rigid 
body displacements, for, although owing to the double integration an arbitrary 
function of a and an arbitrary function of /3 will appear in /iQ, these will be 
determined by (l7), except for functions of a or ^, which make its left-hand side 
vanish identically. The only possible arbitrary terms in /^ Q . are therefore given by 
AQ = a A (cosh a + cos ^) + B(cosh a — cos /3)4-Ca sinh a + Da sin /3, or 

Q == Ar^ + aB 4- CV + Da: 

where r is the distance from the origin. These give rise to terms in ii, v corresponding 
to motions of pure translation and rigid body rotation about the origin. 



§ 3. The Stress-Function. 

Turning now to the consideration of the possible forms for the stress-function in 
these co-ordinates, we note that the differential equation (8) can readily be solved by 
the ordinary method, and that its general solution is 

A^ = e''0i(a+i/3) + 6^\/>2(a-fi/3)4-6'Y3(a---^/3) + 

If we seek a solution of the type hx = f{^) cos /i/3 or f{a) sin n^, (8) shows that 
the differential equation for /(a) is 

|1 2(«^+l)/,+n^-2n^-f l)/(«) = 0, 

the solution of which is 

/ ( a) = A^ cosh (n 4- 1 ) a -f B,j cosh (n — 1 ) ot 4- C,, sinh (n -h 1 ) a -f D,, sinh (n ■— 1 ) a, 
unless 71 = or 1. In the latter case we have 

/ (a) = Aj cosh 2a + Bi -I- Ci sinh 2a -f Dja, 
and when l^ = 

f (a) = Ao cosh a -I- Bya cosh a ~h Co sinh a -I- Doa sinh a, 
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If we now seek solutions for which hx is a multiple of sinh na or cosh na, we find 
the following solutions which are not included above — 

hx = (E cos /3 + F sin ^S+G cosh a+ H sinh a) /3. 

Since any constant multiples of x, y and any constant may be added to x without 
affecting the stresses, it follows from (2) that any multiples of 

sinh a, sin /3 or cosh a —cos /3 

may be added to hx- This allows us to take the coefficients of cosh a, sinh a and sin /3 
as zero. We have then the following general expression for hx : — 

hx = (E cos /3+F sin /34-G cosh a + H sinh a) ^ 
+ (Bo cosh a + Do sinh a) a 
+ ( Ai cosh 2a + Bi + Ci sinh 2a + Dja) cos /3 
-f ( A'l cosh 2a + C'l sinh 2a + D\ol) sin /5 



z' 



00 



+ 2 <1 

n~2 



\. 



[A„ cosh {n+ 1) a + B„ cosh (n— l) a4-C„ sinh (r^+ 1) a 

+ D^ sinh (n—l) a] cos n^ 
+ [A' cosh (n + 1) a + B'„ cosh (rz,— l) a + C'„ sinh (n + l) a 

+ D^^ sinh (n — 1 ) a] sin n/3. 



r^ 



. . . (21) 

We have now to determine whether the displacements corresponding to this stress- 
function are single valued or not. The function (/?.Q) is easily obtained .by simple 
integration from (16), and the arbitrary functions thus appearing can be determined 
by the aid of (17). We have 

- TTlf (^Q) = (E cos /3+F sin ^ + G cosh a + H sinh a) a 

— (Bo cosh a + Do sinh a) /3 

- ( Aj sinh 2a + Cj cosh 2a 4- D'i/3) sin ^ 
+ (A'i sinh 2a + C'l cosh 2a~Di/3) cos /3 

[A'„ sinh {n -f 1) a + B'„ sinh (?i-l) a + C'^ cosh (r?,+ 1) a 

+ D'^ cosh (^ — 1 ) a] cos n^ 

-[A^ sinh (^ + 1) a + B„ sinh (n- 1) a + C„ cosh (n+ 1) a 

+ D„ cosh (7^-- 1) a] sin n/3. 



+ 2 ^ 

n^2 



> 



J 



. . . (22) 

It is clear, from the general expressions for hx and AQ, that the only terms which 
can possibly give rise to many- valued displacements are 

hx = (E cos /3+F sin /3 + G cosh a + H sinh a) ^ 

+ (Bo cosh a + Do sinh a + Di cos /3 + D\ sin ^) a, 

yoL. cox^i. — A, 3 Q 
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and the corresponding terms in /iQ 

-~^AQ = (E cos ^ + F sin ^-f Gcosh a + H sinh a) a 

X + 2./W, 

— (Bo cosh a + Do sinh a + D^ cos /3 + D\ sin ^) /3. 

From (19) and (20) we may now find the corresponding displacements u, v. Each 
of these is found to contain a multiple of the many-valued co-ordinate /3. Equating 
the coefficients of these terms to zero we have the following relations : — 

E + G = 0, Bo+Di = 

mF~-(x + 2m)Do =^ 0, /xH + (x + 2/>t)D^i = O.J ' '-' • 

We shall now show that these early terms correspond to the resultant of the forces 
and couple applied over the boundaries. For this purpose we shall require the following 
elementary forms of the stress-function :— 

(1) For an isolated force X applied at the origin in the direction of the x-axis 

^ = — (27r)"'^X [yd—vx log r) 
where r, 6 as usual denote polar co-ordinates and v = jui/{\-\-2jui). 

(2) For an isolated force Y applied at the origin in the direction of the ^-axis 

X ={27r)~^Y {x6 + vylogr), 

(3) For a point couple of moment L applied at the origin in a positive sense 

^ = --.(27r)-'L0. 

(4) For a centre of pressure radiating uniformly from the origin 

X = log r, 

Inserting the relations (23) necessary to ensure single-valued displacements our 
early terms become 

hx = G (cosh a— cos /3) /3 + /3o (cosh a— cos /3) a 

-h F (^ sin ^ -f VOL sinh a) + H (/3 sinh a — va sin /3) 
or 

-., X = <^Gr/3-l-aBoa-hF(x/3-f-i^|/a)4-H(?//3 — i/Xa) (24) 

Now aG/3 = aG(0i — ^2) and hence this term represents a couple of moment 27raG 
applied at a = 00 and an equal and opposite couple applied at a = — 00 . The term 
a^oOi represents two equal and opposite centres of radial pressure at these same points, 
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We also have 





:?7rA/^^™~4— — - f27rqH 



Il{y/3-'vxa) == 'SL{{y^a)ei—px\ogr^}-'lI{{y-a)e2-vxlogr2} 

This corresponds to a force 27rH applied at a = + oo parallel to the ir-axis and an 
equal and opposite force applied at a = — oo 
(thus forming a couple of moment 47raH) and 
point couples each of moment 27r(xH applied 
at these same points (see fig. 2). 

Finally olsso 

F x^ + vyot.) = F {x9i ~\-v{y + a) log rj 

-F {xe2 + v{y-a)\ogr2} 
— avF log rirg. 

This corresponds to forces each equal to 27rF, p j p 

acting at the points a = ± oo and each directed 

towards the origin, together with two equal like centres of uniform pressure at the 
same points. This brings to light a new solution corresponding to the last term. 
Expressed in our co-ordinates we have 

log rjrg = 2 log (2a) — 2 log (cosh a— cos ^8), 
and the corresponding form of hx is, apart from constants, 

hx = (cosh a— cos ^) log (cosh a— cos /3), 

It is easily seen that this can be expanded in a Fourier series which is included in 
our general expression for Ax, but that the expansion is different on opposite sides of 
the line a = 0. For this reason we shall find it convenient to include a term of this 
form whenever the region under consideration includes parts above and below the 
axis of x, i.e.^ when it is bounded by two circles neither of which encloses the other. 

It will be noted that, taken together, the early terms allow for the most general 
resultant forces acting over the two circular boundaries enclosing the two points 
ot=+oOj a=— CO, subject to the condition that the forces acting over the two 
boundaries considered together form a system in equilibrium. If it is desired to 
investigate problems for which this condition is not satisfied we can readily obtain the 
necessary additional solutions. They will be 

X = {y-\-a)e^-vx\ogr^ 

X=' xe^-^v{y-\-a)\ogr^y , . (25) 

2 Q 2 
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corresponding to forces and couple applied at a = — oo^ and similar terms in Qg, logrg 
corresponding to forces and couple applied at a = +00. The corresponding forms of 
hx ^^^ be expanded in series which are included in our general form, but here again 
the expansions are different on opposite sides of a = and diverge for a = 0, ^ = 
together, i.e,, at infinity. This divergence corresponds to the obvious fact that forces 
or couples must be applied at infinity to maintain equilibrium. 

Owing to difficulties of this kind we shall find it convenient to insert the 
appropriate terms corresponding to the resultant force and couple over a boundary 
and to investigate the stress-function corresponding to the remaining applied forces 
which will be in statical equilibrium for each boundary. 

Let us write for brevity 



(f)^ {a) = A„ cosh (^ + 1) a + B^cosh (n — l) ^4 C^sinh (n + l) a + D,,sinh (^— l) a 
i/'ni^) = A^^cosh {n+l)a + W^cosh. (n— l)a + C'^sinh {n+l)a + 'DJ Binh{n — l)a, 

if n ^ 2 and 

^1 (a) = Ai cosh 2a + Bi + Ci sinh 2a 

\fr^ (a) = A\ cosh 2a + 0^1 sinh 2a. 



-> 



> (26) 



y ...... (27) 



Setting aside the terms corresponding to the resultant forces and couples over the 
separate boundaries we have 

hx = {Boa + K log (cosh a — cos /3)} (cosh a —cos /3) 

+ 2 {^^(a) cosn/3 + '^«(a) sinn/3} . . . . . . . (28) 



71=1 



where the term in K may be omitted when the region considered lies entirely on one 
side of the line 

§4. Boundary Conditions. 

Let us consider a plate bounded by two curves a = a^, ag. We may suppose 
ai > a^ and a^ > 0. Then, if ag > we have a finite plate bounded internally and 
externally by circles which are not concentric, if a^ < we have an infinite plate 
containing two circular holes, and by suitably choosing the values of a, ai, a^ we can 
make the circular boundaries in either case of any desired radius and centre distance. 
In particular if ag = 0, we have a semi-infinite plate bounded by a straight edge 
and containing a circular hole. Suppose that such a plate is in equilibrium under 
given normal and tangential forces applied over the boundaries a = a^, ag, so that we 
are given over a = ai, 



aa/3 = ^0+ 2 (a„ cos n^ + \ sin 7i/3), 

1 

aaa = Co+ 2 (c„ cos n^ + d^ sm n^), 



?" • • • » • . \ZiUJ 
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while over a = otg we have similar expansions in which a^^ a^, &^, Cq, c„, d^ are replaced 
by a 05 <x „5 6 „, Co, c ^, a ^. 

If the tractions applied over the circle a = ai are statically equivalent to forces 
X, Y at its centre, and a couple of moment L, then 



Jo l da daj 

Y = ri^ I- + ^ t^l <^i8, 

Jo I da daJ 



L = --T 



a^ r a/3 d^ 



Jo 



sinh ai J cosh a — cos ^ 

The coefficients of aa, ^^ can readily be expanded in Fourier series. We have, in 
fact, since ai > 0, 

dx sinh a sin B c, Z _««, • o 

da (cosh a — cos p) i 

3m = _^(cosh « cos /3-1) ^ _2^ |, ^^_„, ^^^ 
da (cosh a— cos p) i 

and 

00 . 

sinh ai (cosh ai— cos /3)""^ = 1 + 2 2 e"'*"'' cos 9^/3. 

1 " 

Substituting these and the expansions for aa, a^ in the expressions for X, Y, L, and 
integrating, we have 



00 



X = 27rSm(a„-4)e-''% 

1 

Y = - 2x 2 w (6„ + c„) e-""', 

1 

00 

L = — 27rg cosech^ai 2 a^e~'^'^\ 

The corresponding components of the resultant of the forces applied over a = as 
can be obtained in a similar way. We must, however, remember that in this case 
the forces act from that side of the boundary for which a < ag, whereas in the case 
of the first boundary they acted from the side for which a < aj. We obtain, if a^ > 0, 



00 _ 



00 

X' = -27r2?i(«'„-c?'„)e- 



QO 



Y' = 27rSn(6'„+c'„)e-''^ 



00 

L' = 27ra cosech^ag 2 a'^e""^, 

1 
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If ag < there are some differences of sign owing to the different Fourier expansions 
for the direction cosines. We have 

1 

00 

Y' = -27r2%(&'„-c'„)e»^ 

1 

00 . 

U = —^ira cosechVg 2 a'^e'*"^ 
Hence, if the forces acting on each boundary are statically in equilibrium, we have 



QO CO 

2 n {a„-d^) e-""' = 0, Hn (&„ + c„) e"""' = 0, 



1 

00 



1 



> ■ ■ 



2 a„e-''"> = 0, 



2 a'^e*""'' = 0. 



with, if 01.2 > 0, 



QO 



QO 



2n(a'„-c?'„)e-»'^ = 0, 2w(6'„+Oe-"'^ = 0, 



• • 



or, if otg < 0, 



00 



00 



2 n {a'„ + d'„) e""' = 0, 2 w (6'„-c'„) e""^ = 0, 



• (30) 



• (31) 



(32) 



We will now show that it is possible to determine a stress-function of the form (28) 
which gives the appropriate stresses over a = a^, ot2> ^^^ which gives no stress at 
infinity if the region considered extends so far. 

By the aid of (6) we can calculate the stresses corresponding to the stress-function 
(28). We obtain 



2aaa = K(l— 2 cosh^ a) — 2Bo sinh a cosh a + 29^1 (a) 
+ 2 (K cosh a-f-Bo sinh a) cos /3— K cos 2/3 

' [{n+ 1) (^-1-2) 0^+1 (ot) -2 cosh a (n^-l) (pn M 

+ {n—l){n—2)(p^_i{a)'] G0sn8 
+ [{n+1) {n-i-2) V-„^.i (a) -2 cosh a{n^-l) ^r^ {a) 

-\-[n-l){n-2) x}r^_, (a)] sin n^ 

, —2 sinh a [^'^ (a) cos ^^-}-'\^'^ (a) sin n0]. 



~\ 



+ 2 ^ 



y 



J 



and 



2aoL^ = -^'i (a) — 2 (K sinh a + Bo cosh a) sin ^-{-Bq sin 2^ 



00 



+ 2 



n — \ 



[(% + 1) i/-'„+i (a) -2 cosh an-v!-'„ (a) + (*i- 1) ^'„_i (a)] cos n8, 
- [(ri + 1) /„+] (a) -2 cosh aft^'„ (a) + (w- 1) ^'„_i (a)] sin n^S.' 
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Identifying these with (29)5 we have the following relations from which to obtain 
the coefficients : — 



fpi (ai) = 2cq + 2B0 sinh % cos otj + 2K cosh^ % — K 

1.2.3.^2 (%) = 2ci— 2 (K cosh ai + Bo sinh %) +2 sinh % <f>\ {%) 

2.3.4,^3 {%) —2 cosh % . 1.2.3.^2 (^i) = 4c2 + 2K + 4 sinh ai/gW 

• • * • • * • • « • •'• • •'* • >-• •• • « 

n(n+l)(^4-2)^„+i(aj)-2coshai(n~l)(n)(n4-l)^„(ai)4"(r^--2)(n 

= 2nc„-k-2n sinh % ^^^^ (%) 



>^33) 



l)W^«-i(ai) 



»(n + l) (« + 2) i/'^+i (ai)-2 cosh ai(»-l) («) (« + l) t/r„(ai) + (re-2) (ra-l) (w) V'n-iCai) 

= 2ndn+2n sinh «! i/^'^ («]) (« = l) 
... (34) 

(?i + l)V^'„+i(ai) -2 cosh ai «^/''„(aj) + (w-l)\^'„_i (aj) = 2a„ {n^ l) 



^ • (35) 



2^'2 (ai) —2 cosh aj /j (aj) = — 26i— 2 (K sinh ai + Bo cosh %) 
S^'s (ai) —4 cosh ai ^'g (a,) +^', (aj) = — 26a+Bo 



% 



>■ 



(36) 



• a • 



(w 4- 1) i>'n+i (ai) - 2 cosh % }^/„ (aj) + (n- 1) /„_i (a,) 



2& 



« 



(«' = 3) 



Writing out equations (35), multiplying by e"""' and adding, we have 



n 



(n + l)V''„+i(ai)e-''"'-n>/^'„(«,)e-^"+^''" = 2 S a,e-^' 



or 



« 



(w + l)V/„+i(ai)-w^'n(ai)e~*" = 26"^' 2 a^e 



-paj 



p-0 



. . . (37) 



Now, in yirtne of (30), we may write the right-hand side of this 



00 



00 



2e""' 2 a^e-^"' 



It Zj C^mj.*.6 5 



r = l 



'n+r^ 



and since 2a,, cos n^ is supposed convergent this tends to zero as n increases. Hence 
from (37) we see that the limit of V^^^+i (a,)/^^„ (a^) as n increases is e■-"^ and hence 
the functions ^^^(otj) are finite for all Takes of n and tend to zero as n increases, if 
the resultant couple acting on a = cti Tanishes, 
Multiplying (37) by e^*"^ and adding, we have 



'n—l q 



g=op=o 



'P^ 



2 2 S a„e<'*-^>°', 



p = q-p 
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which, on effecting the summation with regard to q^ leads at once to 

^V^^wC^-i) = 2 cosech cti 2 a^sinh (n— J9)ai . . . , . . (38) 

for n^l. 

Treating (34) in a similar way we have 

B(ti+l)(fi+2)i/.„^i(%)--(^--l)(n)(fi+l)i/-^(%)e^«^ 

= 26""' 2j9[d;p+smhai^/.'p(ai)]e-^'. . . (39) 
We can readily show from (35) that 

% n 

sinh «! S pe-P^^xlr'^ia,) = ^n (n+l) {^'„+i(ai)-e~V«(ai)}e"""'- S pa^e-^S 

and hence (39) may be written 

n{n+l) (w+ 2) x/.„+i {ai)-{n~ l) {n) (w + 1) ^„ (a,) e""' =. n (»+ l).[>/^'„+i (aj) + e-"'V-'„ (ai)] 

+ 26""" 2 J?(dy-a^)e-^^ 

p = 1 

As in the case of ^'^,^(^1)5 we can show that the right-hand side tends to zero as n 
increases if conditions (30) are fulfilled. Hence y^„ (^i) is finite for all values of n and 
tends to zero as n increases, and we have 

n~l q 

{fh-1) (n) {n+ 1) f, (oti) 6^^-^^^^ = 2 2 2 p {d^-a^) e^^^-^^"^ 

+ :E q (g+ 1) e*«' [V^'.+i (aO-e-^'V^'g (^i)], 
which, on reduction, leads to 



n {n^ — l) \f^^ {ai) = 2cosechai 2 {(s— jD)apCosh(n— p)cx] 



+ (p(l^— a^ coth ai) sinh (fi— p) %} . , (40) 

for n^2. Equations (34) do not determine -i/^j (ai). 
From (36) we have 

2^'2 (%) 6-"^-^^i (oti) e""^^^ = ^'1 (ai)-26-~"^ (K sinh aj + Bo cosh %)~26l6'-"^ . (41) 
and II fh — — ^ . 

[n + l) /„+j (a^e-""'- V» («i) e-<"+""' = "A'l (ai)-2Ke- sinh a^-B^ 






and hence, if the sequence <p^ (aj) is to converge for large values of n, we must have 

^^j (oti) = Bp+2Ke"""^ sinh aj-f 2 2 6^e"^"^ ..,,., (48) 

i> = 1 
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From (41) and (42) we have, if n ^ 2, 



?i— 1 
n(^^,(ai)e<«-^^«' - ^^Ja,)~Bo+ Z (^\ (ai)-2K6--» sink ai--B,) e'^^^ 



n—1 q 

-2 2 2 &j,e <'*-">'"■, 

Q- = 1 p = 1 

from which we obtain 



n sinh a^ ^'„ (a^) = (^\ (aj)— Bq) sinh naj — 2K sinh (ti— l) % sinh oti 



— 22 6p sinh (n — p) % (44) 



© = 1 



for n^2y while <{>\ (a^) is given by (43). 

Finally we have from (33), omitting the first equation of the series, i£ n ^ 2, 

n {n + l) {n + 2) e~^^^<j>,^, {a,)-{n^l) (n) (n + l) e-^-^^>^^<^, [a,) 

= -2 (Bo + K) e-"" sinh a^ 



+ 2 2 p (c^+sinh a^ /^ (otj)) 6~^''\ 

jO = 1 



By the aid of (36) we can reduce the right-hand side to 

n 

from which it appears that ^^^i (ai) is finite and tends to zero as n increases, provided 
that the resultant of the applied forces over a = oci is zero. We then obtain for 

^ = 2, 

n [n^—l) sinh a^ ^^ (a^) = {(p\ (ai) — Bo) {n cosh ?^al— coth aj sinh noti} 

— K {{n—l) sinh 77,ai — (n + l) sinh {n--2) a^} 



n~l 



+ 2 2 {(pCp+&^ coth oLj) sinh(n— jp)aj 

— {n—p)bp cosh {n—p)oLi}, . . , . (45) 

while 

2^i(ai) == 2Co + Bo sinh 2ai + K(2 coshV-1) (46) 

It appears that equations (38), (40), (43), (44), (45) and (46) give the values of 
i>n («i)? ^n (^i)? i>^ny («i)> i^^n (^i) ^OY 71 ^ 1 iu torms of Bo, K and the given coefficients 
a^y &c., with the exception of yfr^ (a^). 

Now we have only assumed % > in order to establish the convergence of these 
ftinctions, and hence the corresponding functions of ag will be given by the same 
formulae with a'^, V^, d^, dJ^ substituted for %, 6„, c^, d^, provided that the conditions 
for convergence are satisfied. It may be shown that the new conditions of con- 
vergence are identical with (30) and (31), or (30) and (32), according as a^ > or < 0, 

YOL. CCXXI. — A, 2 B 
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The formula for (j>^ [a^) given in (43), which is itself a condition of convergence, will, 
however, be replaced by 

00 

^'i(a2) = Bo + 2Ke''^sinha2+2 2 &V% (47) 

-?! = 1 

if 062 < 0. 

From (26) we see that the coefficients A^, B^, 0„, D^ for n^2 are determined from 

(pniot-i), (pn{^2)'. '(j>^n{^i)^ <p\{^2), ^^^ slmikrly A.^^, B'^, C'„, D'^ are determined from 

The values of ^i(ai), ^1(^2), ^\(ai), (l>i[oi^ will give four equations to determine the 
three constants A^, Bj, C^, and the condition that they shall be consistent gives one 
relation between Bq and K. The values of ^'^(ai), \lr\{a^ determine the two constants 
A^, C^i, and ^^1(^1), -^1(0^2) are not otherwise determined. 

We have thus just sufficient equations to determine the coefficients in (28) with 
the exception of Bq, K, between which we have found one relation. If a2>0, so that 
the region considered lies entirely on one side of the axis a = 0, we may take K = 0. 
If on the other hand a2<0 the condition that the stress shall vanish at infinity, which 
is hx-^0 when a, /3->0, gives one more relation between the coefficients, so that in 
either case Bq, K are determined. 

We may therefore adopt the following method :■ — Insert terms of the type (24) or 
(25) corresponding to the resultant force and couple on each boundary, and calculate 
the residual stresses over the boundaries. These will now form systems in statical 
equilibrium over each boundary, and we have shown how to determine an appropriate 
function of the form (28). 

The problem of finding the appropriate stress-function for given tractions over the 
boundaries might have been approached by investigating the values of hx and its 
normal gradient on the boundaries, on the lines developed by Michell.^ The direct 
method which we have adopted is, however, in most cases simpler in our particular 
co-ordinates. 

There is an exception to this rule, namely, when a boundary is free from stress. In 
this case the boundary conditions assume a very simple form. From (6) we have 



and 



7^ (hx) = const = />, say (48) 

(cosh a — cos /3) — ^ {hx) — ^i^ ^^^-{hx) + (i08h. a (hx) = p sinh a 

op op 



the solution of which is readily found to give 

hx = p tanh a + o- (cosh a cos ^— 1) + t sin ^ . . . . . (49) 

on the boundary considered. 

■^ 'Proc. London Mathematical Society,' vol. xxi,, 1900, p. 100, 
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The relations (48) and (49) are the necessary and sufficient conditions that a 
boundary a = constant should be free from stress. The constants p, or, r are 
Michell's three constants of the boundary. 

§ 5. A Cylinder or Pipe with Eccentric Bore. 

In this section we will consider the problem of a cylinder, whose cross-section is 
bounded by two non- concentric circles, which is subject to a uniform normal pressure 
over its internal surface and a different uniform normal pressure over Its external 
surface. By Filon's theorem of generalised plane stress precisely the same analysis 
will give the average stresses in a plate of the same section under the same applied 
forces. 

Let the boundaries of the cross-section be defined by a = ai for the internal 
boundary and a = ol^ for the external boundary. Then a^, cl^ ^1*^ positive and ai>a2. 

Let the applied pressures be Pi, Pg respectively, so that aa =— P] on a = a^, aa = — Pg 

on a = otg and a/3 = on both boundaries. 
Let us assume 

IfiX = BoOt (cosh a— cos/3) + (Ai,cosh 2a + Bi-|-Ci sinh 2a) cos /3. 

Calculating aa, ^^8, by means of (6) and applying the boundary conditions, we find 
the following values for the constants : — 

Bo = 2a M (P1-P2) cosh (aj-ag) 

Ai =-aM (P1-P2) sinh (ai + aa) 

Ci = aM (P1-P2) cosh (ai + ag) 

Bi = aM {Pi cosh (aj — ag) sinh 2a2— P2 cosh (aj — otg) sinh 2ai + (Pi + P2) sinh (ai — as)} 

where, for brevity, we have written 

M = |- cosech (a^— ag) {sinh^ai + sinh^ag}""^ 

The most important aspect of the problem is the value of the stress /3/3 in the boundaries, 
for it is upon this that the strength of the cylinder will depend. This is most readily 
determined by (7), and we find without difficulty 

aa— ^i8 = 4M(Pi — P2)(cosh a — cos /3) {sinh (ai + ag— 2a) cos ^— sinh a cosh (a^ -as)} 
so that on a = ai 

)8y8 =— P^ + 4 (Pi— P2)M (cosh ai— cos /3){ sinh (ai — ag) cos/3 + sinh aiCosh (ai — ag)} (50) 
and on a = a2 

Js^ = — Pg- 4 (Pi — P2) M (cosh ag— cos /3){sinh [a^ — oL^ cos 8— sinh cl^ cosh (ai — ai)} (51) 

2 R 2 
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In order to investigate these results further we will consider separately the cases 
when the cylinder is subject to either internal or external pressures. There is no 
greater difficulty in the consideration of the general case, should the necessity arise, 
except that the formulae are correspondingly longer. 



A Cylinder under Internal Pressure. 
If we put P2 = 0, we have on the external surface 

^^ =— 4P1M (cosh Ota— cos^){sinh (ai — ots) cos 8— sinh otg cosh (ai — ag)}, 

if c?i, d2 denote the distances of the circles oti, ay from the origin, ri, rg their radii and 
d the distance apart of their centres, so that d = d2—di we may show from (2) that 

di = a coth oti, c?2 = a coth otg 

r^ = a cosech oti, rg = a cosech otg 

d, = {ri^r,'^d')lU, d, = {r,'^r,'-^d')/2d 

a'=:{r,'^{r,-\-dY}{ri^{r,^dY}lU\ 

By means of these relations we can reduce the expression for /3/3 to the form 

— _ 2P,nH^/(n-2c? cos 0)'-{r,'-dJ} 

From this and the obvious inequality d! <r2-ri we easily see that — 

(1) The numerically greatest stress is when /3 = tt, Le., on the line of centres at the 

thinnest part of the cylinder. This is always a tension if Pi is positive and is 
given by 

(n^ + r/)(r/-n^-2r,(^ + c?^) ' ' * ^ ^ 

(2) If the centre distance is greater than half the external radius there is minimum 

stress at the points corresponding to cos ^ = rjid. This is always negative 
when Pi is positive and we have maximum compressions equal to 

^V^rHn'-^d'Y (53) 

(ri^+r/){r/-(ri + c?)^}{r/-(ri-c?)^} . . v ; 

This is always numerically less than the maximum tension. There is a 
secondary maximum at ^ = 0, Le,, on the line of centres at the thickest part of 
the cylinder, which is equal to 

2Pi(r/ + n'-2r,d!--cZ^) /g^x 
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(3) If the centre distance is less than half the external radius, we have, in addition 
to the maximum tension (52), a minimum at ^ = given by (54). There are 
no other maxima or minima and the stress decreases steadily from its value at 
the thinnest part of the cylinder to its value at the thickest part. 

On the internal surface we have 

^/3 = — Pi4-4PiM (cosh ai— cos ^) {sinh (oti — ag) cos /3 + sinh ai cosh (oti — otg)}, 
or, expressed in terms of the radii and centre distance, 

-_ p 2fir/{(r/-(^y-n^(n + 2(^cosm .„. 

Hence it may be shown that 

(1) If the centre distance is greater than one-half the internal radius the maximum 

stress in the internal surface occurs at the points corresponding to 
cos 13 = — rj/2c? and is 

' (n'+r/){r/-(n-(^r} {r,'-{n+df} ' • ' • • ^''''^ 

(2) If the centre distance is less than one-half the internal radius the maximum 

stress is at /3 = tt, i.e., on the line of centres at the thinnest part of the cylinder. 
It is 

(3) The minimum stress is at /3 = 0, the point where the line of centres meets the 

internal boundary at the thickest part of the cylinder. It is 

__p + 2'P, r,'{r2' + n'+2r,d-d') .^gx 

^ (ri^+r/)(r2^-r/4-2ri(i-<i^) ^ ^ 

This may be shown to be essentially positive if P is positive so that, as would be 
expected, the internal boundary is everywhere in a state of tension. 

A Cylinder under External Pressure. 
Putting Pi = in (50) and (51) we have on the internal surface 

/3j8 =— 4P2M (cosh ai— cos/3) {sinh (ai — ag) cos /S+sinh a^ cosh (ai — ag)} 

2P,r/ {(r/^cg^)^-r/ {r, + 2d cos ^f} . . 

(n^+r/){r/-(n-d!)^}{r/-(n + d)^} ^ ' 

and on the external surface 

/3^ == — P2+4P2M (cosh ag— cos /3) {sinh (oti — ag) cos /8— sinh ag cosh (oti— otg)} 

p 2IV,^ {r/ {r,-2d cos jSY-jr-^-d r} /..x 

- -^2 (r{'+r,'){r,'-{r,-dy}{r'-{r, + df} ^ ' 
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Hence if the centre distance is less than half the internal radius the compression in 
the inner surface decreases steadily from a maximum at the thinnest part of the 
cylinder to a minimum at the thickest part ; otherwise there is a minimum at each of 
the points and maxima at the points corresponding to cos ^ = —r^jid. Similarly if the 
centre distance is less than half the external radius the compression in the outer 
surface decreases steadily from a maximum at the thinnest part of the cylinder to a 
minimum at the thickest part ; if the centre distance exceeds this value the 
compression is a maximum at each of these points and minima at the points corre- 
sponding to cos ^ = r^l'Zd. 

If in these results we put d = 0, we have, for a concentric tube under internal 
pressure, tensions at the inner and outer surfaces which are respectively 

while for a tube under external pressure the compressions at the inner and outer 
surfaces are respectively 

2r/ -p r/+r/ 



^ p 

These are the well-known formulae for thick tubes. 

§ 6. A Semi-infinite Plate with a Cibcular Hole Subject to a Uniform 

Normal Pressure. 

If in the results of the last section we put ag = and Pg = 0, we have the solution 
for a semi-infinite plate containing a circular hole, which is subject to a uniform 
normal pressure, and bounded by a straight edge which is free from stress. 

We have on the boundary of the hole 

^/3 = --P1 + 2P1 cosecy oLi (cosh^ a^— cos^ /3) 

and on the straight edge 

^/3 = — 2P1 cosech^ oti (1 — cos /3) cos ^. 

If r is the radius of the hole, d the perpendicular distance of its centre from the 
straight edge, and x the distance measured along the straight edge from the foot of 
the perpendicular, 

d = a coth a, r = a cosech ai, d^—r^ = a\ 
and 

X ~ a sin /3/(l — cos 0). 

We have therefore on the straight edge 

Tb- ~4P r'{x'--d'-VT') . . 
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This has a maximum tension at the symmetrical point {x = O) of magnitude 

A-py/id'-r'), . . . . . . . . . . (62) 

At the points x = ±\/{(P—r^) it vanishes, and then becomes a compression which 
reaches a maximum value at points at distances ±^3{d^—r^) on either side of the 
foot of the perpendicular from the centre of the 
hole, which is numerically equal to one-eighth of 
maximum tension. 

The stress round i)he circular hole may be 
represented by a simple geometrical construction. 
If in fig. 3 the centre of the circular hole is C, Q 
is any point on the circle, and CA the perpen- 
dicular drawn from C to the straight edge, and if 
Q denote the angle QAC, we easily see that 

tan ^ = sin ^ cosech otj, 

and the stress round the circular hole is 




/3^ = P,(l4-2tan2 0) (63) 



Fig. 3. 



Hence the stress is the same at points Q, Q^ which lie on the same ray through A. 
The stress is minimum at the points nearest to and most remote from the straight 
edge, where it is a tension P numerically equal to the applied pressure. Thus at 
these points the stress is the same as it would be in the absence of the straight 
boundary if the plate were infinite. The maximum stresses are at the points of 
contact of the tangents drawn from A the circular boundary. At these points its 
value is 



Pi 



d^—r^ 



(64) 



The maximum tension in the circular boundary is equal to the maximum tension in 
the straight edge if d = v^Sr. In this case each is equal to 2Pi. If the distance of 
the hole from the straight edge is greater than this value the maximum tension is at 
a point on the circular boundary ; and if it is less, the maximum stress tension is at 
the symmetrical point on the straight edge. On the other hand, the point of 
maximum difference of principal stresses is on the straight edge or the circular 
boundary, according as d is greater or less than \/2r. This suggests a simple 
method of determining whether, for a particular material, rupture occurs at the point 
of greatest tension or at the point of greatest stress-difference. If a circular hole is 
bored near the straight edge of a uniform plate, so that the distance of its centre 
from the edge is greater than ^2 and less than -v/3 of the radius, and a uniform 
radial pressure is exerted over the hole in any convenient way and increased until 
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rupture occurs, the crack will begin on the straight edge according to the greatest 
tension theory, and on the edge of the hole if the greatest stress-difference theory- 
holds. 

It will be noted that the stresses produced will become large if the hole is near to 
the straight edge. The formulae are so simple that it is hardly worth tabulating 
their numerical values, but a single example will serve as an illustration. If the 
shortest distance from the hole to the straight edge is one-tenth of the radius of the 
hole, the maximum tension in the straight edge is 19*5 times the pressure in the 
hole. 

§7. A Semi-infinite Plate Containing an IJnsteessbd Circular Hole and 
Under a Uniform Tension Parallel to its Straight Edge. 

Let the circular boundary be defined by a = oti, so that if r is its radius and d the 
distance of its centre from the straight edge, 

T = a cosech a^, d = a coth ai, d/r = cosh aj. 

At a distance from the hole the stress-function may be taken as x = ¥^y^ where 
T is the tension, so that, if a > 0, 

h-^Q = ^aT sinh^a/(cosh a —cos /3) 

= iaT sinh a U + 2 i ^-'^'^ cos ?^^) . . ... . (65) 

We have to add to this a stress- function which gives no stress at infinity and no 
stress over a = 0, and is such that the complete stress-function gives no stress over 

We may omit the term in K in (28), since in this case the region considered lies 
entirely on one side of a = 0, and clearly the required stress-function is even in /3. 

It may readily be seen that the condition that aa and a^ shall vanish over a = is 
satisfied by (28) if cp^ (O) = and (p^ (o) = for ^ ~ 1, and hence from (26) and (27) 
A^ + B^ = and (r2, + l) C^^H- (n+l) D^ = 0. We may therefore take for our complete 
stress- function 



hx = aT 



|- sinh all -^2 2) e~"^ cos n/3\ + BQa (cosh a -cos /3) + Ai (cosh 2a- 1) cos ^ 

I n = 1 



GO 



+ 2 <; 

11 = 2 



An [cosh (r?, + 1 ) a — cosh {n— 1 ) a] ^ 

y cos n/3 
+ E,;[(r2,— l) sinh {n-{-l) a — {n-{-l) sinh {n—l) a] 



. . (66) 



At infinity a = 0, /3 = the first series diverges, but may of course be replaced by 
the alternative form in (65). If the second series converges it is clear that at infinity 

X = Xo- 
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We must now choose the coefficients in (66) so as to satisfy (48) and (49), and 
there is no difficulty in finding the following values for the coefficients : — 

Ai = |-e-^"'sech 2ai, Bo = sech 2ai (67) 

A — — n^ Bmh^oLi—n sinh % cosh ai + e"""'"' sinh nai rr.r,\ 

'"' 2 {m^h%al—n^ sinh^a^} 

-p, _ n sinh^oti 



2 {sinh^noti— n^ sinh^oti} ' 



Substituting in (66) we have for the complete stress-function, 



Ax = «T 



a sech 2ai (cosh a — cos 0)-\-^ sinh a + sech 2ai cosh (2a] —a) sinh a cos ^ 
'n sinh otj sinh (a— %) sinh na ] 

^ . —sinh a sinh n (a— a^) sinh n%' cos n^ 



n = 2 



sinh%ai— n^ sinh^oti 



.(69) 



We may now calculate the stress ^^ in the boundaries by means of (6). We find, 
on the circular boundary a = oti, 



00 



^^1 = 2T (cosh oti— cos /3)-|smh a^ sech 2ai+ 2 M^ cos n/3^ . . : (70) 

I 71=2 J 

where 

M — ^(^--l) Binh{n-\'l) ai—n{n-\-l) sinh(n-~l) oti /^.v 

2 {sinh^noti— n^ sinh^ai} 

The stress in the straight boundary cannot be directly determined from (69), for it 
is found that the resulting series diverges for a = 0. We can, however, find without 
difficulty from (66) that when ct = 0, 

^^^ = TJl + (l-cos^)iP«cosn^| ...... (72) 

where P. = 4nA.. ^ ^ ' ^ 

The series in (70) converges only slowly, unless a^ is large, and for convenience in 
computation we may transform it by separating the more slowly converging part. 

Let 

M,, = 2n(nsinh oti— coshoti) 6~"'*''' + N^, ...... (73) 

and we readily obtain 

2 (cosh ai-cos 0) ^ n{n sinh oti- cosh a^ e"""^^ cos nB ^ 1 - . 2Bmh^%sm ^ 

n = i (cosh oti — cos /3)^ 

Substituting in (70) we have 

iS r=i 2T [l ^ sinh^ a^ sin^ jB 
\ (cosh ai — cos 0f 

-f 2T (cosh ai-cos^)-! sinh ai sech 2aj + 2e~^'*^ cos /3+ 2 N^cosn^}^ . (74) 
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If 6 is the angle between the radius to the point a^, ^ and the perpendicular to the 
straight edge, then 



sin 6 



sinh a, sin ^ 
cosh oti — cos/3 



and if ctj is large (74) reduces to ^^i = T(l + 2 cos 29), which agrees with the known 
result for a hole in an infinite plate and gives compression numerically equal to T at 
the extremities of the diameter parallel to the tension, and tensions equal to 3T at 
the extremities of the perpendicular diameter. 

The numerical values of the coefficients P^, N^ are given in Tables I. and II. 
respectively. It will be noted from Table I. (that, as ai increases, Pg tends to become 
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large compared to the other coefficients. Hence, when the hole is at a considerable 
distance from the straight edge, the stress in the straight edge approximates to 

T{l"-C cos 2/3(1 -cos /3)} 

where G is a small positive constant. 

This shows that the stress in the straight edge is a minimum at the mid-point, 
increases to a maximum as we move outwards, then diminishes to a second minimum, 
and finally increases steadily to the value T at infinity, where /3 = 0. 

In fig. 4 we have plotted the graphs of the stresses in the boundaries for a case in 




t&O" 



Stress ill straight edge. ------ Stress in circular boundary. 

Fig. 4. 

which the hole is fairly near to the straight edge, a^ ™ 0'8, for which the shortest dis- 
tance between the two boundaries is approximately one-third of the radius of the circle. 
It will be noted that the general character of the stresses is not afl:ected by the proximity 
of the straight edge. It will be remembered that when the hole is at a great distance 
from the straight edge there are maximum stresses of 3T at the extremities of the 
diameter perpendicular to the straight edge, with points of maximum compression 
numerically equal to T lying between. For a^ ~ 0*8 we find that the maxima occur 
at the same places but are increased, the increase being more marked at the point 
nearest to the straight edge, where the tension is 4*366T, while its value at the point 
most remote from the straight edge is 3'266T. The stress in the straight edge also 
maintains the same general character as it exhibits when the hole is at a great distance 
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from the straight edge. Here again, however, the maxima and minima are 
accentuated. The minimum at the central point has decreased and has become a 
compression numerically equal to 1*956T. 

It appears that for the range of values vvhich we have investigated the maximum 
stress is on the circular boundary at the point of nearest approach to the straight edge. 
Its value for different values of the ratio of the distance of the centre of the hole from 
the straight edge to the radius of the hole, together with the stresses at the centre of 
the straight edge and at the most remote point of the circular boundary, is shown in 
Table III. 
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It will be noted that when the hole is very near to the straight edge, so that the 
two boundaries are separated only by a narrow connecting piece, the stress in this 
piece consists of a very large tension on the inside and a numerically slightly less com- 
pression on the outside. Hence, as might be expected from general considerations, 
the stress in this narrow connecting piece is a bending moment accompanied by a 
certain amount of tension. 

These results may be compared with some experimental results recently obtained by 
Prof CoKER and Messrs, K. C. Chakko and Y. Satake^^ by optical means. These 
deal with the stresses in a strip of finite width under tension with a circular hole 
centrally placed, whereas we have considered the case of a semi-infinite plate with a 
circular hole near its straight edge. The problems are therefore not quite comparable ; 
but as in each case the critical region will clearly be near the minimum section between 
the hole and a straight boundary, the two problems may be expected to exhibit the 
same general characteristics. For the strip of iinite width it is found that there is 
maximum stress in the circular boundary at the points of nearest approach to the 



* ' Transactions of the Institution of Engineers and Shipbuilders in Scotland,' vol. Ixiii., Part I., p. 33, 
1919. 
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straight edges and minimum stress at the points of the straight edge immediately 
opposite the centre of the hole. Moreover as the radius of the hole is increased in 
proportion to the distance of its centre from the edges of the strip these maxima and 
minima become more pronounced. In all the cases examined experimentally the 
minimum stress in the straight edge remains a tension, but Prof. Coker surmises that 
if the radius of the hole were still further increased in proportion to the width of the 
strip this minimum stress would become a compression. All these results agree 
qualitatively with the theoretical results established in this paper for the semi-infinite 
plate, and allowing for the difference in the two problems they may be taken as a 
substantial experimental verification. 



